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It is well known that the charge current in a conductor is proportional to the applied electric 
field. This famous relation, known as Ohm's law, is the result of relaxation of the current 
due to charge carriers undergoing collisions, predominantly with impurities and lattice vi- 
brations in the material. The field of spintronics 1 5 , where the spin of the electron is ma- 
nipulated rather than its charge, has recently also led to interest in spin currents. Contrary 
to charge currents, these spin currents can be subject to strong relaxation due to collisions 
between different spin species, a phenomenon known as spin drag 6 } This effect has been 
observed for electrons in semi-conductors 7 and for cold fermionic atoms 8 9 , where in both 
cases it is reduced at low temperatures due to the fermionic nature of the particles. Here, we 
perform a transport experiment using ultra-cold bosonic atoms and observe spin drag for 
bosons for the first time. By lowering the temperature we find that spin drag for bosons is 
enhanced in the quantum regime due to Bose stimulation, which is in agreement with recent 
theoretical predictions^. Our work on bosonic transport shows that this field may be as rich 
as transport in solid-state physics and may lead to the development of advanced devices in 
atomtronics^ll 



1 Body 



Ohm's law can be written as j = E/p, showing the simple linear relation between the applied 
electric field E and the current density j. The electrical resistivity p comprises a wealth of infor- 
mation on the material and its careful measurement over the last century in solid-state physics has 
introduced new intriguing phenomena, like high-temperature superconductivity 12 and the integer-^ 
and fractional 14 quantum Hall effect. For the case of spin currents in two-component cold-atom 
systems with pseudo-spin up and down a similar relation applies^ 

Js = A (1) 

Ps 

where the applied force F\ on particles with spin up, mass m and density n leads to a drift velocity 
difference Av between the two spin species. This results in a spin-current density j s = nAv, 
which is inversely proportional to the spin resistivity p s . Because of the absence of impurities and 
an ionic lattice the spin resistivity p s = my /n is for ultra-cold atoms completely determined by 
the spin-drag rate 7. This rate is a measure for the transfer of momentum via interactions between 
the two spin species. The similarity between Ohm's law and Eq. ([T]) can be understood from the 
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Figure 1 : Schematical representation of the experimental methods - a) Collision of two atoms 
with different spinP. The mass current j m in this process as shown by the horizontal arrow is 
conserved, whereas the spin current j s as shown by the vertical arrow is not conserved. In this 
example the spin current is reversed, b) Schematic representation of the measurement: Two spin 
species (blue, yellow) are prepared and a force acts on only one spin species (yellow). Due to the 
spin drag the other spin species (blue) is dragged along. 

analogy between their underlying physical mechanism. In both cases the external force or field 
leads to an acceleration, which on average is interrupted after a time r and this leads to a finite 
current. For spin currents in ultra-cold atoms the particles with one spin collide with particles of 
another spin and these collisions do not conserve the spin current, as shown in Fig. [T^. This is 
the microscopic origin of spin drag. It is important to emphasize that in electronic systems spin- 
drag effects also exist, but are usually small compared to the effects of impurities and phonons^B, 
whereas for ultra-cold atoms they are the only effect and therefore dominant. 

To observe spin drag in trapped ultra-cold atomic quantum gases, we set atoms with pseudo- 
spin up into motion with respect to atoms with pseudo-spin down and observe, how atoms with 
spin down are dragged along with the atoms with spin up. In our first method, the constant force 
method, we apply a constant force Ff on only the spin species up as shown in Fig. [IJ5. After a 
short period a constant drift velocity difference At? has developed between the two spin species, 
where the ratio of the applied force F\ and Av is proportional to the drag rate 7. This method is 
well suited for a large drag rates, but for weak drag a steady state may not be achieved before the 
two clouds are spatially separated. Therefore for small drag rates we apply our second method, the 
oscillation method. In this method the centers of the two spin species are first separated from each 
other and then start to oscillate in the trap. The decay time of the relative oscillation is proportional 
to the drag rate. For both methods we have to detect the change in relative position of the clouds 
as a function of time. We apply a Stern-Gerlach technique^ in a direction perpendicular to the 
applied force to spatially separate the spin species and image the clouds with absorption imaging. 

In the experiment we load up to 4.6 x 10 8 23 Na atoms into a cigar-shaped optical far-off- 
resonant trap (FORT) with trap frequencies of o; ra d/27r = 835 Hz in the radial and u; ax /27r = 3.5 
Hz in the axial direction. The trap is tight in the radial direction to obtain a large density, thereby 
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producing a collisionally opaque and hydrodynamic sample in the axial direction. The temperature 
is between 2-8 /iK and is always kept above the critical temperature for Bose-Einstein condensa- 
tion. The atoms in the trap are initially all in one spin state. For stability and to maximize the 
interspecies collisions we then flip about 50% of the atoms to another spin state, creating an equal 
mixture of spin species up and down (See Methods). 

In Fig. [2] we show an example of the constant force method. The force is applied on the spin 
species up such that this cloud starts to move to the left and due to spin drag the spin species down 
is dragged along. Careful observation shows that the upper cloud of spin species up is skewed: 
The left wing of the cloud is longer than the right wing. This is a first, direct indication of spin 
drag, which is shown more clearly in Fig. [2}) in the profiles along the pulling direction through the 
center of each cloud. The clouds are fitted to skewed Gaussian distributions as shown in Fig. [2}d 
to determine the center of mass. As indicated with the two vertical lines in Fig. [2J5, the difference 
in position between the two clouds is much smaller than the width of the clouds such that local 
effects can be neglected. 

To obtain the spin-drag rate this measurement is repeated for different magnitudes of the 
force and different durations during which the force is applied. In all cases we measure the relative 
displacement as shown in Fig. [2}). The result for different duration t is shown in Fig. [3] After a 
short time the difference in position between the two clouds is linear in time, indicating that there 
is a constant drift velocity Av between the two spin species. This is direct evidence for spin drag 
as the drag leads to friction between the two spin states. In the absence of drag the difference in 
position would show a quadratic dependence on time due to the ballistic motion of the atoms. The 
curvature for short times (t < 5 ms) in Fig. [3] is due to the fact that the two clouds are initially not 
in a steady state. However, the clouds are collisionally opaque, i.e., they are in the hydrodynamic al 
regime, and the two clouds reach quickly a steady state. In our analysis we allow for this by 
solving the coupled equations of motion for the centers of the two clouds 16 and the results are used 
to extract 7 from the data. 

Equation [T] is only valid for small applied force Ff, such that the response of the clouds 
is still in the linear regime. Two of us^ showed that for ultra-cold atoms to be in the linear- 
response regime Av should be well below the thermal velocity v t h = /rem, where ks is 
the Boltzmann constant and T the temperature. In the experiments t> t h is of the order of 6 cm/s, 
whereas typical values for Av are of the order of 1 mm/s. Thus we expect that our experiments 
are in the linear regime. In the inset of Fig. [3] we have tested this assumption by measuring Av as 
a function of the acceleration Ff/m. The line is a linear fit that shows good agreement with the 
data, which indicates that this measurement is indeed in the linear regime. Since we conducted this 
measurement with the lowest number of particles and thus lowest drag rate, this proves that all our 
measurements for the constant force method are in the linear-response regime. 

As mentioned, for lower drag rates we applied the oscillation method. Here, only 3-7 x 10 6 
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Figure 2: Measurement for the constant force method - A single measurement for the constant 
force method on a cloud with 1.5 x 10 8 atoms at 4 pK. The spin species up is accelerated with 1.01 
m/s 2 to the left, a) Pseudo color image of the measured absorption along a line-of-sight through 
the two clouds. The upper cloud is in the spin up state, whereas the lower cloud is in the spin down 
state. The separation between the clouds results from the Stern-Gerlach technique after the spin 
drag to facilitate the imaging of the individual clouds. The difference in width of the clouds in 
the vertical direction develops during the detection stage, c) Absorption of the two clouds in the 
horizontal direction, where the lower curve is for the upper cloud (spin up) and the upper curve 
is for the lower cloud (spin down). The two curves are obtained by making a cut through the 2D- 
image of Fig. [2^ at the height of the respective clouds. The red/orange curves are for spin state 
up/down and represent a fit with the skewed Gaussian distributions for the individual clouds. The 
skew in the two distributions is a direct result of the drag of the upper cloud exerted on the lower 
cloud. Notice that the separation between the two clouds is small compared to the width of the 
clouds. 
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Figure 3: Experimental results for the constant force method - Relative displacement Ax be- 
tween the position of the two clouds as a function of duration t that the force is applied. The data 
points are indicated by the red dots and the blue line is a fit to the data using the solution of the 
coupled equations of motion of the two clouds. After a short period (t < 5 ms) Ax is linear in 
time, indicating that there is drag between the two clouds. All data points are offset due to an 
artefact in the Stern-Gerlach technique, which shifts after the spin drag the axial position of the 
two clouds with respect to each other by a constant amount. The inset shows the drift velocity Av 
as a function of the strength of the force Ff. The straight line is a fit to the data, which shows that 
the experiments are in the linear-response regime for the drag, where Eq. ([TJ) holds. 
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atoms are loaded into the FORT. Atoms with spin up are spatially separated from atoms with spin 
down by exerting with a constant force on the spin species up for a period, after which the force 
is switched off. The atoms in spin state up start to oscillate in the trap and drag along the atoms 
in spin state down, which are initially at rest. In Fig. [4] we show the relative displacement between 
the two clouds as a function of time after releasing the atoms. Since the drag rate 7 is much larger 
than the axial trap frequency co ax , the relative oscillation is overdamped and the decay rate is given 
by Pd — 2o; ax 2 /7, which allows us to extract the spin-drag rate 7 from the measurements. Even for 
the lowest drag rates the oscillation remains overdamped. 

As explained earlier, spin drag is caused by collisions between atoms in different spin states. 
It is therefore to be expected that the drag rate 7 should be proportional to the interspecies collision 
rate 7^. For a non-degenerate gas Stringari and Pitaevskii 17 find 7 = 1 7-1-4.. The factor 2/3 < 1, 
independent of temperature, reflects the fact that for the drag not all collisions contribute equally 
because only the velocity component in the direction of the force is relevant. In Fig. [5] we plot 
our results for the drag rate as a function of the interspecies collision rate, which is given for the 
trap by 7^4. = (n)a^v re \. Here (n) is the density of one spin species averaged over the trap, 074 
is the interspecies collisional cross section and v rc \ the relative velocity between the two colliding 
atoms averaged over the velocity distribution. As shown in Fig.[5J our results in the non-degenerate 
regime are in good agreement with the result of Stringari and Pitaevskii 

For larger phase-space densities the gas can no longer be considered as classical and quantum 
statistics becomes important. For our Bose gas collisions are enhanced at low temperatures due 
to Bose stimulation. This should be contrasted with fermionic systems, where the drag rate is 
suppressed at low temperatures due to Pauli blocking^. Collisions are enhanced due to two effects. 
First, the central densities of the gas in a trap with a fixed number of particles and temperature 
increases due to Bose enhancement beyond its classical value. Second, the scattering process is 
enhanced due to Bose stimulation to the final states. In this regime we no longer expect to find a 
temperature-independent proportionality factor between the drag rate and the interspecies collision 
rate. Therefore, in Fig. [6] the drag rate is shown as a function of the fugacity z = exp(/x//csT), 
where /i is the chemical potential of the gas. For 2 <C 1 the gas is in the classical limit, whereas 
for z — 1 the gas starts to Bose-Einstein condense. Van Driel et a/™ showed that the drag rate 
scaled as 7/T 2 is a function of the fugacity z only. In Fig. [6] we plot 7/T 2 as a function of the 
fugacity z, where z and T are determined from the experiment (see Methods). For a classical 
gas 7/T 2 depends linearly on z and the classical result is indicated in the figure with a straight 
line. For small z our measurements agree with this result, but for larger values of z the drag rate 
increases above the classical rate. Evaluating the expression of van Driel et al^ for the drag rate 
of an inhomogeneous sample, we find the Bose enhancement of the drag rate, as indicated in Fig. [6] 
by the blue line. As can be seen from the figure, our results are in good agreement with theory. 
As a measure of the accuracy of the theory we use the reduced chi-squared parameter x 2 ed , which 
becomes unity for a large number of observations when the data agrees with the theory. For the 
classical result we obtain x 2 ed = 6.10, whereas the calculations including Bose enhancement yields 
X 2 ed = 1.58. We emphasize that here Xr C d * s use d to compare our experimental results to a theory 
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Figure 4: Experimental results for the oscillation method - The relative displacement Ax in 
position between the two clouds for the oscillation method as a function of time t after the release of 
the atoms. The relative oscillation is overdamped, in which case the decay rate (3d of the oscillation 
is given by (3d = 2u; ax /7 with u; ax the oscillation frequency in the trap and 7 the drag rate. The 
red points are the data points and the blue curve is a fit to the data assuming exponential decay. 
All data points are offset due to an artefact in the Stern-Gerlach technique, which shifts after the 
oscillation the axial position of the two clouds with respect to each other by a constant amount. 
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Figure 5: Drag rate versus the interspecies collision rate in the classical regime - Drag rate 7 
as a function of the interspecies collision rate 7^. Here the results for the constant force method 
are indicated with black dots and the results of the oscillation method are indicated by green dots. 
The red line represents the prediction of Stringari and PitaevskiP^ and in this range of 7^ agrees 
with our experiment results. 
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without any adjustable parameters. 

It is important to realize that the experiments presented here can only be performed when 
the sample is in the hydrodynamic regime. In this case the number of collisions during one trap 
oscillation is much larger than one. To quantify this we can define a hydrodynamicity parameter 
= 7coi/w ax , w ^ m 7coi the total collision rate including inter- and intraspecies collisions. For the 
current experiments we have reached 7 co i = 1.4 kHz at an axial trap frequency of co a _ x /2ir = 3.5 Hz, 
yielding a hydrodynamicity of > 60. This large hydrodynamicity allows us to quickly reach a 
steady state. This can for instance be seen in Fig. |3} where a steady state is reached after 5 ms 
corresponding to only a few collisions. 

In conclusion, we have studied spin drag and determined the spin-drag rate for a gas of ultra- 
cold bosonic atoms. We find that in the quantum regime spin drag is Bose enhanced and up to 
almost a factor of two larger than the classical value, which is in good agreement with a recent 
theoretical prediction 10 . This shows that the measurement of spin drag is a strong precursor of 
Bose-Einstein condensation. Our results pave the way for studies of the transport properties of 
degenerate bosons which are very different from fermionic systems, such as electrons in ordinary 
conductors, and give a complementary picture of transport. Our results are among the first that 
study transport properties of ultra-cold atoms and in the long term are important for the develop- 
ment of atomtronic devices that are based on transport of ultra-cold atoms^R 

2 Methods Summary 

Sodium atoms are trapped and cooled in the \F = l,m = — 1) = | f) state (pseudo-spin up) in a 
magnetic trap, before transferred to a FORT. In the FORT a radio-frequency sweep is applied, such 
that half of the atoms are transferred to the state \F = 1, m = 0) = | \) state (pseudo-spin down). 
The atoms in the FORT are tightly confined in the radial direction and the sample is hydrodynamic 
in the axial direction. The force on the atoms is applied using a magnetic field gradient, which 
only acts on the spin up atoms. After the experiment the two clouds are spatially separated using 
a radially inhomogeneous magnetic field and subsequently imaged using absorption imaging. The 
spin-drag rate is extracted by solving coupled equations for the position of the individual spin 
species and extracting the relative displacement Ax from the solutions where the spin-drag rate is 
a free parameter. This solution can be compared to the experimental data to extract the spin-drag 
rate. 

The inter- and intraspecies collision rates are determined using the measured number of 
atoms in each cloud and the temperature. Combined with the known scattering cross section o and 
the trap frequencies in the radial and axial direction the inter-, intraspecies and total scattering rate 
can be determined. 
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Figure 6: Drag rate versus the fugacity in the degenerate regime - Scaled drag rate 7/T 2 as a 
function of the fugacity z = exp(/i/fcs). The dots represent the data points, whereas the curves 
show the results for the spin-drag rate using the classical results of Stringari and Pitaevskif^ (red 
curve) and the quantum results of van Driel et alP^ (blue curve). For small z the gas behaves 
classically and the data agrees with both models. For large z (z > 0.5) the data deviates from the 
classical model and is in agreement with the quantum model. The results show that the measure- 
ment of the spin-drag rate is a precursor for the phase transition of Bose-Einstein condensation at 
z= 1. 
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3 Methods 



Preparation The samples were prepared in a magnetic trap by standard laser and evaporative 
cooling methods of sodium atoms in the \F = l,m = —1) = | f) state (pseudo-spin up). There 
are two major improvements with respect to other experiments concerning the particle numbers: 
Spin polarizing atoms before the transfer from the magneto-optical trap to the magnetic trap and 
decompressing the trap axially at the end of the evaporative cooling stage 19 . The final trap frequen- 
cies are axially Lo ax /2n = 1.5 Hz and radially aj rad /27r = 116 Hz with up to 8 x 10 8 atoms. After 
evaporatively cooling the atoms to a few /xK the FORT is adiabatically turned on in 400 ms with a 
trap depth of about 50 /jK, an axial trap frequency of u slx /2tt = 3.5 Hz and a radial trap frequency 
of u! Ta d/2n = 835 Hz, leading to 4.6 x 10 8 atoms in the trap. The lifetime is more than 7 s. The 
atoms in the \F — 1, m = 0) = | I) state (pseudo-spin down) are produced with a radio-frequency 
sweep from below the bottom of the magnetic trap: 2.6 MHz to 2.71 MHz in 30 to 50 ms. The 
fraction of atoms in the \F = 1, m—+l) state is estimated to be less than 10% and the same for all 
measurements and does not influence the drag rate. 



Experiment and Detection The force on the spin species up is produced by a magnetic field 
gradient using anti-Helmholtz coils with a constant offset magnetic field such that the zero crossing 
of the field is not in the area of interest. The detection of different spin species is achieved by 
turning on the radial magnetic confinement for 5 ms only and then ramping down the FORT in 
2 ms for radial cooling. The radial cooling is performed to assure spatial separation of the two 
spin species. After this time the spin species down falls freely due to gravitation while the spin 
species up is still radially trapped. An additional time of flight of 5-10 ms before imaging is used 
depending on the density and temperature of the clouds. The imaging is performed by sending a 
probe beam through the two clouds and detecting the absorption on a CCD camera. 



Analysis The axial position difference is determined by fitting the clouds with a Gaussian distri- 
bution multiplied by 1 — Erf (cr(x — x )), where Erf is the error function, x the central position in 
the Gaussian distribution and a a measure of skewness. From the fits the number of atoms of each 
species and the temperature are deduced. The latter is deduced from the length of the clouds when 
the force is applied for 1 ms only. Each measurement is repeated three times for error reduction. 
Temperature and total number of atoms are also measured separately by detecting a non-perturbed 
cloud expanding freely in 30 ms time of flight and fitted with a Gaussian distribution. Especially 
for large clouds this is necessary since the anharmonicity of the trap in the axial direction causes a 
systematic error in the temperature measurement. 



The differential equations describing the system 



are 



1 dV(x+) F t 

x t --fAx + = -4 (2) 

m dx^ m 

il + lAi + l d ZM = o ( 3) 

m axi 
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with the axial center-of-mass position of spin species up and down, respectively, Ax = x^—x^, 
m the atomic mass, V(x) the axial trapping potential and F-j- the applied force on the spin species 
up. Here drag rate 7 = T(Ax) / (/t A m n tot ^in) * s me spin-drag rate with T(v) the drag force 
on the two clouds, n tot the total density of the atoms, and the fraction of atoms with spin up 
and down, respectively. Neglecting the effect of imbalance in T(v), this expression corrects to 
leading order for non-equal spin mixtures and neglects higher-order corrections. These equations 
are numerically integrated using the initial conditions x^(0) = x±(0) = and from the results the 
difference in position Ax as a function of time t is determined. It can be shown by using similar 
equations for three spin states that the drag rate between two of these species as defined above is 
not influenced by the presence of a third species, if the three drag rates are the same. 



The classical interspecies collision rate is obtained by integrating the interspecies collision 
rate over the trap and we fincP^ 

7t4 = WWrel = ttH^T ' ( } 

where v ie \ is the relative velocity between the two atoms, N is the total number of atoms and 
er-1-4, = 47ra^ 2 the interspecies cross section. For 23 Na the scattering length is = 2.80 nm^. 



In the oscillation method the separation between the clouds is induced using the same mag- 
netic field as in the constant force method and the displacement of spin up is well within the range 
where the FORT is harmonic. The solution of the equation of motion is 

Ax(t) = Ax(0) exp ( J (5) 



with Ax(0) the initial position difference and f3 d the decay rate. From /3 d we can deduce the drag 
rate 7 using 

2 ^ax 

7 = ~f^- (6) 

PA 

To determine 7 theoretically we evaluate the microscopic expression for the spin-drag rate includ- 
ing Bose enhancement given by van Driel et aP. 
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